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Abstract: 

We consider a Higgs boson coupled to gluons via the five-dimensional effective operator 
H\xG iW G^ v produced by considering the heavy top quark limit of the one- loop coupling 
of Higgs and gluons in the standard model. We treat H as the real part of a complex 
field 4> that couples to the selfdual gluon field strengths and compute the one-loop 
corrections to amplitudes involving eft, two colour adjacent negative helicity gluons and an 
arbitrary number of positive helicity gluons - the so-called i^-MHV amplitudes. We use 
four-dimensional unitarity to construct the cut-containing contributions and the recently 
developed recursion relations to obtain the rational contribution for an arbitrary number 
of external gluons. We solve the recursion relations and give explicit results for up to four 
external gluons. These amplitudes are relevant for Higgs plus jet production via gluon 
fusion in the limit where the top quark mass is large compared to all other scales in the 
problem. 
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1. Introduction 



As the time for physics at CERN's Large Hadron Collider (LHC) approaches there is a great 
need for precision calculations of Standard Model processes. The eagerly anticipated new 
physics signals typically result in complex multiparticle final states that are contaminated 
by Standard Model contributions. Isolating the signal therefore relies on precise theoretical 
calculations of background rates. 

At present, the rates for such processes can be estimated relatively easily at lead- 
ing order (LO) in parton level perturbation theory. A wide range of tools are available, 
for example ALPGEN [1, 2], the COMPHEP package [3, 4], HELEC/PHEGAS [5, 6], 
MADGRAPH [7, 8] and SHERPA/AMEGIC++ [9, 10] which compute the amplitudes 
numerically 1 and provide a suitable phase space over which they can be integrated. How- 
ever, the predicted event rates suffer large uncertainties due to the choice of the unphysical 
renormalisation and factorisation scales so that the calculated rate is only an "order of mag- 
nitude" estimate. In addition, there is a rather poor mismatch between the "single parton 
becomes a jet" approach used in LO perturbation theory and the complicated multi-hadron 
jet observed in experiment. 

While these problems cannot be entirely solved within perturbation theory, the situ- 
ation can be improved by calculating the strong next-to-leading order (NLO) corrections. 
The uncertainty on the rate may reduced to the 10%-30% level by including the next-to- 
leading order (NLO) corrections, which may themselves be large when new channels are 
accessed. In addition, the additional parton radiated into the final state allows a better 
modelling of the inter- and intra-jet energy flow as well as identifying regions where large 
logarithms must be resummed. 

The ingredients necessary for computing the NLO correction to an particle process 
are well known. First, one needs the tree-level contribution for the n + 1 particle process 
where an additional parton is radiated. Second, one needs the one-loop n particle matrix 
elements. Both terms are infrared (and usually ultraviolet) divergent and must be carefully 
combined to yield an infrared and ultraviolet finite NLO prediction. 

The real emission contribution is relatively well under control and can easily be auto- 
mated [1-4, 7-9, 11]. The infrared singularities that occur when a parton is soft (or when 
two partons become collinear) can then be removed using well established (dimensional 
regularisation) techniques so that the "subtracted" matrix element is finite and can be 
evaluated in 4-dimensions [15]. 

The bottleneck in deriving NLO corrections for multiparticle processes is computing the 
one-loop amplitudes. The standard approach is to compute the relevant Feynman diagrams 
using a variety (or combination) of numerical and algebraic techniques. Much progress has 
been made in this way, numerical evaluations of processes with up to six particles have 

^OMPHEP, MADGRAPH and SHERPA/AMEGIC++ evaluate sums of Feynman diagrams. However, 
both ALPGEN and HELAC/PHEGAS are based on off-shell recursion relations [11-14]. 
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been performed [16-18]. However, one always observes large cancellations between the 
contributions of different Feynman diagrams and the result is generally far more compact 
than would naively be expected. This is a strong hint that more direct and efficient ways 
of performing the calculation exist. The simplicity can be realised using on-shell methods, 
where the cancellations due to gauge invariance and momentum conservation are already 
present, to compute the amplitude. 

In a series of pioneering papers, Bern et al [19, 20] developed the use of on-shell methods 
at loop level by sewing together four-dimensional tree-level amplitudes and using unitar- 
ity to reconstruct the (poly)logarithmic cut constructible part of the amplitude. A key 
feature of the unitarity method is that the tree-level helicity amplitudes are often very sim- 
ple. However, the rational terms, that are produced when computing non-supersymmetric 
amplitudes in D = 4 — 2e dimensions, were difficult to obtain. 

Recently, on-shell methods have received renewed attention with Witten's proposal 
of a duality between J\f = 4 supersymmetric Yang-Mills and a topological string theory 
[21]. New analytical methods, the MHV rules [22] and on-shell recursion relations [23], 
have been derived and have been extremely successful in computing tree-level helicity 
amplitudes. Both of these techniques can be proved using simple complex analysis and a 
knowledge of the factorisation properties of the amplitudes on multi-particle poles [24, 25]. 
Although developed for in the context of multigluon amplitudes, they have had wide ranging 
applications including processes with massive coloured scalars [26, 27], fermions [28-31], 
Higgs bosons [32, 33] and vector bosons [28, 34]. A more complete set of references can be 
found in refs. [35, 36]. 

Some of these on-shell ideas have been applied to one-loop amplitudes including the 
application of maximally-helicity- violating (MHV) vertices [37], a more generalised unitar- 
ity [38] using complex momenta and the use of the holomorphic anomaly [39]. Together 
with previous results [19, 20], these new methods have led to the complete analytic expres- 
sions for the cut-containing contributions [40-42] for one-loop QCD six-gluon amplitude. 
Very recently an efficient method for direct extraction of the integral coefficients exploiting 
complex momenta has also been proposed [43]. 

However, these techniques are also intrinsically 4-dimensional and, although providing 
effective ways of computing the cut-constructible (poly) logarithmic part of the amplitude, 
suffer from the same limitations as the older unitarity based methods and miss the rational 
part. These "missing" rational terms can be obtained directly from the unitarity approach 
by taking the cut loop momentum to be in D = 4 — 2e dimensions [44] and there have been 
recent developments in this direction [45-48]. 

On the other hand, the rational part is essentially tree-level-like in containing only poles 
in the complex plane. One can therefore attempt to isolate these terms using an on-shell 
recursion relation in an analogous way to tree level amplitudes. This is the "unitarity on- 
shell bootstrap" technique which combines unitarity with on-shell recursion [49-51]. Other 
Feynman diagram based methods have also been developed for calculating the rational 



- 2 - 



terms directly [52-56]. Altogether, it has already been possible to calculate the full set 
of complete 6 gluon helicity amplitudes as well as closed forms for n-point amplitudes for 
specific helicities [57-59], and the complete set of six photon amplitudes [60, 61]. 

As mentioned earlier, the MHV rules have also been successfully applied to tree level 
QCD amplitudes involving a massive colourless scalar, the Higgs boson, in the large top 
mass limit [32, 33]. This is achieved using a decomposition of the Higgs field into its selfdual 
and anti-selfdual components, <f> and (ft . The purpose of this paper is to apply the new on- 
shell methods to one- loop calculations of the Higgs plus n-gluon amplitudes. Amplitudes 
of this kind have been considered in the context of on-shell recursion in relations by Berger, 
Del Duca and Dixon [62] for the finite helicity configurations that vanish at tree level, the 

0- "all plus" and 0- "almost all plus" amplitudes. The cut-constructible parts of the infrared 
divergent amplitudes involving the 4> and an arbitrary number of negative helicity gluons 
were studied in ref. [63]. Here we consider the amplitudes with two adjacent negative 
helicity gluons and any number of positive helicity gluons, the c/>-MHV amplitudes. These 
represent one of the ingredients needed for the complete set of helicity amplitudes for Higgs 
to four gluon process which has been computed using numerical methods in reference [64]. 

Following [51], we will split the calculation into two parts, evaluating the "pure" 4- 
dimensional cut-constructible C n and rational R n parts of the leading colour contribution 
to the one-loop n-point amplitude separately: 

A^ = C n + R n . (1.1) 

For simplicity, we have dropped the leading colour subscript, = A n l \ C n contains 
all of the terms originating in box, triangle bubble loop integrals which are related to 
cut-containing logarithms and dilogarithms (or 7r 2 ) and which are cut-constructible in 4- 
dimensions. In this paper we choose to compute the cut-containing contribution using the 

1- loop MHV rules developed by Brandhuber, Spence and Travaglini [37]. In addition, C n 
contains fake singularities that come from tensor loop integrals. To explicitly remove these 
spurious singularities, it is convenient to introduce a cut-completing rational term, so that 
the "full" cut-constructible term C n is given by 

C n = C n + CRn, (1.2) 

with the corresponding modification to the rational part, 

Rn = Rn — CRn. (1-3) 

Because the rational part contains only simple poles, the aim is to construct this recursively 
using the multiparticle factorisation properties of amplitudes. This means constructing 
a direct recursive term, i?^, by summing over products of lower point tree and one-loop 
amplitudes. By construction, R® encodes the complete residues on the physical poles. The 
cut-completion contribution CR n may also give a contribution in the physical channels 
which would then lead to double counting. These potential unwanted contributions are 
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removed by the overlap term, so that 

£ n = < + 0„, (1.4) 

and the full amplitude is given by 

4 1 ) = C n + CR n + R% + O n . (1.5) 

Our paper is organised as follows. In Sec. |2| we introduce the complex scalar field 
<f> and review the effective interaction that couples it directly to gluons. The relationship 
between eft (and (jr) and Higgs amplitudes is spelt out, along with out conventions for colour 
ordering amplitudes and a selection explicit results for tree-level (^-amplitudes. We address 
the computation of the cut-constructible C n and cut-completion CR n parts of the one-loop 
amplitude in section |3[ The rational contributions are considered in section |I| where we 
establish the on-shell recursion relation and give expressions for the overlap terms O n . 
We solve the arbitrary multiplicity results given in sections || and || for the special case 
n = 4 in section || and give an explicit analytic expression for the A± amplitude. Section || 
is dedicated to a series of checks of our results. We show that our results have the correct 
infrared pole structure, satisfy the correct collinear limits and we study the limit where 
the momentum of the Higgs boson becomes soft. Finally, our findings are summarised in 
section Two appendices are enclosed that define our spinor notation and list the relevant 
one- loop functions that appear in C n . 



2. The Higgs Model 



In the Standard Model the Higgs boson couples to gluons through a fermion loop. The 
dominant contribution is from the top quark. For large mt, the top quark can be integrated 
out leading to the effective interaction [65, 66], 

C 



c 



int 



(2.1) 



In the Standard Model, the strength of the interaction C has been calculated up to order 
0(oii) in [67]. However, for our purposes we need it only up to order 0(a 2 s ) [68], 



C 



6irv 
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1 + — — + ... 

4 7T 



(2.2) 



with v = 246 GeV. This approximation works very well under the condition that the 
kinematic scales involved are smaller than twice the top quark mass Mt. [69-71] 

The MHV structure of the Higgs-plus-gluons amplitudes is best elucidated [32] by 
considering H to be the real part of a complex field <f> = \(H + iA), so that 

C 



£int 
H,A 



HtiG uu + iAtiG w/ *GT 



C 



tr Gsd uv G^ D + 4>HrG 



ASD flU ^ASD 



(2.3) 
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where the purely selfdual (SD) and purely anti-selfdual (ASD) gluon field strengths are 
given by 

Gl v D = \(G» V + *GT) , GZd = \{G^ ~ *G flu ) , 

with 

The important observation of [32] was that, due to selfduality, the amplitudes for <j> plus 
n gluons, and those for (ft plus n gluons, each have a simpler structure than the gluonic 
amplitudes for either H or A. Amplitudes can be constructed for plus n gluons and for 
(ft plus n gluons separately. 

The decomposition of the HGG and AGG vertices into the self-dual and the anti-self- 



dual terms eq. (2.3) means that the Higgs and pseudoscalar Higgs amplitudes are obtained 
from cj) and (ft amplitudes. 

At ] (H, g x \...,gt) = A^ (<f>, <?i \ • • • , 9* n ) + A^ {(ft , g x 1 , . . . , g^ ) , (2.4) 
4 m) (A, g x \ . . . , g x " ) = (</>, 5i 1 , ■ ■ ■ , gt ) ~ A^ ((ft , g x 1 , . . . , g x " ) . (2.5) 

However, parity further relates (f> and (ft amplitudes, 

At\ft,g*\- ■ -,g x n) = (AtH<f>,gI Xl ,- ■ ■ ,9n A ' 1 ))* • (2-6) 

From now on, we will only consider <j) amplitudes, knowing that all others can be obtained 
using eqs. (O)-pJ). 



2.1 Colour ordering 

The tree level amplitudes for a 4> and n gluons can be decomposed into colour ordered 
amplitudes as [72, 73], 

A^(^{k l ,Xi,a i }) = iCg n ~ 2 £ tr(TMD . . . f'W) A^^, <x(r\ ..,**»)). (2.7) 

creS n /Z n 

Here S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj labels the 
momentum pj and helicity Xj of the j th gluon, which carries the adjoint representation 
index otj. The T ai are fundamental representation SU(iV c ) color matrices, normalized so 
that Tv(T a T b ) = 5 ab . The strong coupling constant is a s = g 2 /(4ir). 

Tree-level amplitudes with a single quark-antiquark pair can be decomposed into 
colour-ordered amplitudes as follows, 

An(<f>, {Pi, A;, ai}, {Pj,\j,ij}) (2.8) 
= iCg n ~ 2 J2 (™ 2) ■ ■■T a " (n - 1) )hi n A n ((f>, l\a(2 x \. . . , (n - l)^" 1 ),^) . 

CTGSn-2 

where S n -2 is the set of permutations of (n — 2) gluons. Quarks are characterised with 
fundamental colour label ij and helicity Aj for j = 1, n. By current conservation, the quark 
and antiquark helicities are related such that Ai = — A ra = A where A = ±~. 
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The 1-loop gluonic amplitudes which are the main subject of this paper follow the same 
colour ordering as the pure QCD amplitudes [19] and can be decomposed as [62, 63], 

[n/2]+l 

A${<j>,{k i ,\ 1 a i }) = iCg n E G n . c (a)AU(<t>,a(l^,...,n x «)) (2.9) 

where 

G n; i(l) = Ntv(T ai ...T an ) (2.10) 
G n;c (l) = tr(T ai • • • T^- 1 ) tr(T ac • • • T a ") , c> 2. (2.11) 

The sub-leading terms can be computed by summing over various permutations of the 
leading colour amplitudes [19]. 

2.2 Tree level <p amplitudes 

As noted in [32] the all-plus and almost all-plus 4> amplitudes vanish, 

AW(cj>,g+,g+,g+,...,g+) = 0, (2.12) 
A^(cj ) ,gi,g+,g+,...,g+)=0, (2.13) 

for all n. 

The tree (^-amplitudes, with precisely two negative helicities are the first non-vanishing 
(f> amplitudes. These amplitudes are the </>-MHV amplitudes and general factorisation 
properties now imply that they have to be extremely simple [32] . For the case when legs q 
and p have negative helicity, they are given by 

t 

(12) (2 3)---(n-l, n) (n 1) 
In fact, the expression eq. (|2.14 ) for the 0-MHV n-gluon amplitude has precisely the same 



4 O) (0, 2 2 + , • • • , 9p , ■ ■ ■ , 9^, ■ ■ ■ , 9t) = /inWno v i Zw-TV . ( 2 - 14 ) 



form as the MHV n-gluon amplitudes in pure QCD [74]. The only difference is that the 
total momentum carried by gluons, Pi + P2 + ■ ■ ■ + Pn = —Ptf, is the momentum carried by 
the </>-field and is non-zero. 

There are two 0-MHV amplitudes involving a quark pair, 

_ _ _, (r l) 3 (rn) 

A n {(j),q 1 ,...,g r ,...,q^) = ——— — — - 2.15 

(1 2) (2 3) • • • (n — 1, n) (n 1) 

(v l\ (v 71^^ 

A M , q t,..., 9 -,. = (12>(23) \..' ( L' li „ >( „ 1> . (2.18) 

We note in passing that the tree (^-amplitude with all negative helicity gluons, the 
^-all-minus amplitude, also has a simple structure [32], 

4 o) (0;i-,---,n-) 

4 

= (-11" (2 17) 

1 ' [12][23].->-l,n][nl] ' 1 ' 
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Amplitudes with fewer (but more than two) negative helicities have been computed 
with Feynman diagrams (up to 4 partons) in Ref. [72] and using MHV rules and on-shell 
recursion relations in Refs. [32, 33]. 

The main goal of this paper is the construction of the one-loop </>-MHV amplitude 
with two adjacent negative helicities. For definiteness, we focus on the specific helicity 
configuration (1~, 2~, 3 + , • • • ,n + ). 

3. The cut-constructible contributions 

In a landmark paper, Brandhuber, Spence and Travaglini [37] showed that it is possible to 
calculate one-loop MHV amplitudes in M = 4 using MHV rules. The calculation has many 
similarities to the unitarity based approach of Refs. [19, 20], the main difference being that 
the MHV rules reproduce the cut-constructible parts of the amplitude directly, without 
having to worry about double counting. This is the method that we wish to employ here. 

The four-dimensional cut-constructible part of one-loop amplitudes can be constructed 
by joining two on-shell vertices by two scalar propagators, both of which need to be con- 
tinued off-shell. A generic diagram is shown in figure |] and the full amplitude will be a 
sum over all possible permutations and helicity configurations. In the BST approach the 
propagators are continued off-shell as for the tree level MHV rules and can be written [37] , 

Li = li + Zir]. (3.1) 

Loop integration over the Lj are related to phase space integration over and trivial 
integrals over z\. Note that this approach actually performs all integrals and therefore there 
is a one-to-one identification of cut diagrams and cut-constructible functions. However, it 
relies on being able to perform the phase space integral. For the cases considered in this 
paper, the phase space integrals are known and are directly related to loop functions. 

A generic diagram can be written: 




Figure 1: A generic one loop MHV diagram or unitarity cut. 
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where A-l{R) are the amplitudes for the left(right) vertices and P is the sum of momenta 
incoming to the right hand amplitude. The important step is the evaluation of this ex- 
pression is to re-write the integration measure as an integral over the on-shell degrees of 
freedom and a separate integral over the complex variable z [37]: 

Li^ Li 2 Z\ Z2 

= {z ^tz + ^/ hd%SM(ll)SM(ll), (3.3) 

where z = Z\ — Z2 and z' = Z\ + z 2 . The integrand can only depend on z, z' through the 
momentum conserving delta function, 

$(*)( Ll -L2-P) = <5 (4) (/i -l 2 -P + zr)) = 6&(h -h- P), (3.4) 

where P = P — zr\. This means that the integral over z' can be performed so that, 

V = ^Sf^ /t/ d4z i^^( + )(Z?)^ + )a|)<5W(Z 1 -Z 2 -P) J 4 z a 1 ,-P,-Z 2 )^(Z 2) P,-Z 1 ) 
= (Ai) 2 2m /y/ dL1 ? si4) (-h,l2,P)A L (h,-P,-l 2 )A R (l 2 ,P,-h), (3.5) 



where, 



dLlPS^(-h,l 2 ,P) = ^d%d%5^\ll)^ + \ll)5^{h -l 2 -P) (3.6) 

(2vr) 4 



The phase space integral is regulated using dimensional regularisation. Tensor integrals 
arising from the product of tree amplitudes can be reduced to scalar integrals either by using 
spinor algebra or standard Passarino-Veltman reduction. The remaining scalar integrals 
have been evaluated previously by van Neerven [44]. 



3.1 Pure cut contributions 



The pure cut contribution is constructed by connecting two tree-level vertices and the seven 
independent topologies are shown in Figure |2| Note that the last four topologies helicity 
configurations allow both fermionic and gluonic contributions. Note also that all fermion 
loops always appear in association with a factor of Np, the number of fermion species, and 
a factor of —1. 

Let us consider diagram |2|(a) to begin with. We can take the momenta to be labelled 
from 1 to n around the right hand amplitude so that we consider a cut in the s\ tn channel, 
i.e. when the momentum flowing across the cut is pi + p 2 + • • • + p n corresponding to 
an invariant mass si jn = (pi + p 2 + . . . + p n ) 2 - 2 Other diagrams with this topology are 
accessible by permuting the arguments of the si jn channel. The product of the two vertices 

2 Note that when j < i, Sij = (p» + . . . + p n + pi + ■ ■ -Pj) 2 - 
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Figure 2: The MHV loop diagrams contributing to the <p — > g 1 g 2 g£ . . . g+ amplitude, 
can be written: 

At Au = {12) " 

{hh) 2 {nh){hl 2 ){l 2 l)T\ n a -^{aa + l) 

Applying a Schouten identity to the numerator and using momentum conservation in the 
form h = l 2 + Pi, n w e find, 

A A _ 4(0) ( N(P ljn , Pl ,p n ) Pl.n'Pn Pl.n ; Pi \ ,„ „x 

^ (fc-p^fe+pO* (iT^)2 + ) > (3.8) 

where N(P,p 1 ,p 2 ) = P 2 (pi ■ p 2 ) - 2(P ■ p{)(P ■ p 2 ). This is now written in terms of scalar 
integrals so we can directly use the results of van Neerven [44] to perform the phase space 
integration: 



J d D UPS(-h,l 2 ,P) — 



N(P, Pl ,p 2 ) 
+ Pi) 2 {h+P2) 2 



CF 2isin(7T6) / u 2e |P 2 r e 2 F 1 ( 1, -e; 1 - e; j!/ 2 ^ ) (:.!.!)) 



(4vr)2 e 2 v 1 1 ^ V ' ' 'N(P, Pl ,p 2 ) 

J d D UPS(-hJ 2 ,P)^^ = -^- 2 2is^e)^\P 2 r (3-10) 

J d D LIPS(-h,l 2 ,P) =- -—^L _ — 2zsin(7re)/x 2e |P 2 |- £ (3.11) 
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where the factor cr is given by, 

_ ( ^ni±gizi). 

The final integration is over the z variable. However, the only dependence on z appears 
through the quantity Pi >n 3 so it is convenient to make a change of variables, 

* = «?L (3.13) 

Z p2 _ p2 

to produce a dispersion integral that will re-construct the parts of the cut-constructible 
amplitude proportional to (si )n ) -< S 

/ ^^^ sin ( 7re )i p2 i" £ = 2 ™(- p2 y- ( 3 - 14 ) 

The final result for this diagram then reads: 

= £E A ,o, / _£Y ( Fl L , _ ^ ) + ^ . (3 . 15) 
e 2 V- S i,n/ V V N(P hn , Pl ,p n )J J 



The other "gluon-only" channels (Figs. |2|(b) and (c)) reduce to scalar integrals in the 
same way and we merely quote the results, 



V 2 \ I rn ( 1 1 Pi' PnSl, n -l 

" e; 



+ 2 Fi fl,- e; l-e; 



£ 2 V- S l,n-1/ \ V' ' iV(Pl, n -l,Pl,Pn) 



,n— 1 > Pn j Pn— lj 

~ 2^1 f 1, - £ ; 1 " e; „fi' Pw - iai -- 1 + l) (3-16) 

v u = * A m ( jiLV ( 2Fl (i, _ e; i _ e; . 

e 2 \-si,zJ \ V ^(Pi,i,Pi 5 Pi+i, 
+ 2 Fi l,-e;l-e; 



Fi l,-e;l-e; 



N{Pi^,p n ,pi) 

Pi ■ PiSl,j 

N(P 1:hPl , Pi ) 



- 2 Fi l,-e;l-e;— — . (3.17) 

V N(P 1} i,p n ,p i+ i)J J 

The arguments of these expressions can be straightforwardly permuted to produce results 
for channels ||(a),(b) and (c) for all gluon configurations. 

When considering the channels with alternating helicity configurations around the 
loop we find that even after the Schouten identities have been applied, we are still left with 



^Through a suitable choice of r], one can always ensure that N(P,pi,p2) is independent of z [37] 
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(i 



JVpX / 2tr_(2Pnl) 3 
N ) I 3 S ? 2 (2F.p„)3 



tr (2Pnl) 2 \ . /-, 
si 2 (2P.p n y) + 4 I 1 



JV F \ tr_(2Pnl) 
47V J Si 2 (2P.p n ) 




Figure 3: Decomposition of the MHV diagram of fig. 2(d) contributing to the P2,n channel 



tensor integrals which must be further reduced to scalar integrals by expanding in terms of 
all possible tensor structures. This feature has also been seen in the context of finding the 
cut-constructible part of pure QCD amplitudes and was also addressed applying Passarino- 
Veltman reduction [42]. For diagram ||(d), the s 2jn channel, the presence of tensor integrals 
and fermion loops results in new structures of order 1/e. The result for this diagram is, 



V 



2,n c r 



1+2^1 



+ 2 Fx I 1, — e; 1 — e; 
- 2 Fi ( l,-e;l-e; 



P2 • PnS2,n 
N(P2, n ,P2,Pn) 

Pi ■ PnS2,n 

N(P 2 

Pi ■ P2S2,n 
N(P 2 , n , Pl ,P2) 



1 



N F 
~N 



2tr_(2P 3 ,n-ml) 3 tr_(P 3 ,n-ml) 2 



3sf 2 (2P-n) ; 



i _ JV^ ftr_(J$», n _inl} 



4/V,/ \ si 2 (2P-n) 
where we have introduced the shorthand notation, 

tr-(abcd) = (ab) [be] (cd) [da 



1 - 2e 



1 - 2e 
(3.18) 

(3.19) 



This is better illustrated in figure || which shows the cuts of each integral function that 
appear. Figure Q shows the decomposition of the s 2j j channels (figure |2|(e)) which follows 
exactly the same steps as the previous case. Diagrams ^(f) and ||(g) are analogous to 
diagrams §(d) and ^(e) and can be found be permuting the arguments: 1,2, ...,n — > 
2,1, n, ... ,3. 

When summing over all the possible diagrams we find that the bubble integrals always 
appear in the combination, 

Bub(s) - Bub(i) = £>(e°). (3.20) 

Hence all of the 1/e poles coming from bubble functions vanish leaving the expected combi- 
nation of boxes and triangles proprotional to the tree amplitude[75, 76]. The combination 
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O + i) 



(1 
(1 



/ + i\Y i^iy 



JVpX / 2tr_(2P(z+l)l) 3 _ tr_(2P( ? +l)l) 2 \ , W-i _ Np\ tr_(2P(i+l)l) 
N ) I 3 S ? 2 (2P.p l+1 )3 s\ 2 {2P. Pl+1 y j AN) s 12 {2P.p i+1 ) 



Np\ ( 2tr_(2Pzl) 3 _ tr_(2Pzl) 2 \ . . /-, Np\ tr_(2Pil) 
iV ) \3s 3 12 (2P. Pi f sl 2 (2P. Pi y) ^^\ L AN) s 12 (2P. Pi ) 



Figure 4: Decomposition of the MHV diagram of fig. 2(e) contributing to the P 2 ^ channel 



of bubble integrals can be written in terms of a basis of pure logarithms, 

log(a/t) 



L k (s,t) 



(s-t) 



k ' 



(3.21) 



These logarithmic contributions are not proportional to the tree amplitude, but are multi- 
plied by new spinor structures written in terms of traces. The full, unrenormalised result 
for this specific MHV helicity configuration for n > 3 is thus: 



C n (<M-,2-,3 + ...,n+) = c r 4 0) (<M-,2-,3 + ,...,n+) 



^ (F3 m (s i)ri+i _ 2 ) - F3 m (s i)ri+ i_i)) 



n n+i— 1 



+ 



1ft lb ^ L _L I I 

i=l j =i+3 i=l 

iY P tr_(2P 3 ,i-iil) 3 



_(lP iin (i-l)2) J 



i=4 



^3(^-1,1,5^1) + 



'12 



'12 



N P ti.{lP i)n {i-l)2f iVptr_(2P3,^iil) 2 r , 

"Tj 2 ^2(^-1^,5^1) g - L 2(S2,i; s 2,i-l, 



'12 



'12 



JVptr_(lP i)n (»-l)2) s , iV P tr_(2P 3 ,i-iil) r , > 



+ 



6 S12 
A) tr_(lP <|W (t-l)2) 

N Si2 



6 S12 
/3b tr_(2P 3 ,i-iil) 



A/ 



S12 



£l(s2,i,S 2 ,j-l) 



, (3.22) 



where the one-mass triangle Fg m and box functions F4 are defined in Appendix |B|. For 
convenience, we have introduced 



Pa 



UN - 2N F 



N P = 2(1 



N F 

~~n 



(3.23) 
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Note also that summations of the form Y^ a are understood to vanish when b < a. If we add 
this term together with its complex conjugate with the appropriate momentum re-labeling 
then we find agreement with the known Higgs MHV amplitude in the case of n = 3 [77] . 



Eq. ( 3.22j ) can be rewritten in a form which is more convenient when computing the 
completion and overlap terms, namely 



C n (<M-,2-,3+...,n+) = 
cr4 0) (^l",2-,3 + ,...,n+) 



^ (F3 m (Si,n+i-2) - F3 m (s i) „ +i _i)) 



i=l 



n n+i—1 



n/ u / j F 4 [SiJ, Si + lJ-i; SiJ + l, Si + lj) — - 2_^F 4 (Si,j+2; St,i+1) s i+l 



i+2) 



i=l j=i+3 



i=l 



+ 



n« =2 (aa + l) 



N P 

~6~ 



(lP 4 , n (t - 1)2)(1(» - 1)P, 2 2> - 1)P, 2 2) - (lP i>n (i - 1)2) 



xL 3 (si-i,i,a f) i) + ^ (12) 2 (ip, n (i - l)2)L a (si_i,i,a f) i; 
+^(li J Fb >i _i2)(lPi, i _ii2)((lPi, i _ii2) - (lzP 3li -i2; 



X-M s 2,i, s 2 ,i-i) + (1 2) 2 (lzP 3) j_i2)L 1 (s 2 ,i, s 2) i-iy 



(3.24) 



3.2 Cut-completion terms 



The functions Lk(s,t) that appear in eqs. djj ) and are bubble contributions pro- 



duced by the reduction of tensor box and triangle integrals. They contain unphysical 
singularities as s — > t so it is useful to redefine the cut-containing contribution in terms of 
a new basis which has good behaviour in the various limits. This is achieved at the cost of 
adding some rational terms, 



L l {s 1 t) = L 1 {s,t) (3.25) 
Hs,t)=Us,t) + ^(] + l) (3.26) 

^^^^^(J + J). (3-27) 



The new functions are now free from spurious singularities. Replacing the functions 
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(a) (b) (c) 

Figure 5: Factorisation of one-loop amplitudes in the limit Pf — * 0. Contribution (c) does not 
appear for MHV amplitudes. 



in (3.24) with the corresponding functions we are left with the "completed cut term": 
C7? n (<M~,2-,3 + ,...,r^ 



i+ „+\ - CrNp 



{lPi,n(i ~ 1)2) (1(» - l)Pi, 2 2) ((l(t - l)P i)2 2) - (lPi, n (i - 1)2)) 

1 ( \ 1 



i=4 



X 



+<lii^,i_i2)<lP 1)i _ii2)f<lP 1)< _i*2) - (liP 3 ,i-i2; 



i n i 

+ 



(S2,i - •S2,i-l) 2 \ s 2,i S 2 , 



' t-1 



(3.28) 



4. The rational contributions 

By definition, the rational terms only contain poles in the invariants and therefore are 
amenable to the same type of analysis as used for tree-amplitudes. We will therefore 
calculate the remaining rational terms in AjP((f), 1~ ,2~, ... ,n + ) using on-shell recursion 
relations that have been successful in both QCD [49-51, 57-59] and for the finite amplitudes 
involving the (j) field [62]. 

This approach relies on both the factorisation properties of one-loop amplitudes on 
physical poles and the introduction of a complex shift parameter z to study the behaviour 
of the amplitude in the complex plane. 

We recall the multi-particle factorisation properties of one-loop amplitudes [78] , 

4 1) -^^ + % 2 ^+^%J^ as i?- 0> (4.1) 

it i 

where the subscripts L and R denote the amplitudes with fewer external particles on the 
left and right of the factorising physical channel as shown in fig. || Note that T only 
contributes in multi-particle channels if the tree amplitude contains a pole in that channel. 
MHV amplitudes do not have multi-particle poles and hence this term is absent. 
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Following the discussion in the earlier sections, we divide the amplitude into its cut- 
constructible and rational terms 

A^=C n + R n . (4.2) 



Applying eq. (4.1) in the rational sector we find (for MHV amplitudes), 



i i 

However, as it stands, R n contains unphysical poles which cancel with unphysical poles 



in C n . Therefore, we add the rational function CR n (the cut completion term of section 3.2 ) 
to C n and subtract it from R n . CRn is chosen such that it cancels the unphysical poles in 
C n (and thus, simultaneously in R n ). This means that we should try to set up a recursion 
on the physical poles with the expression R n = R n — CR n , 

To develop the recursion, a convenient analytic continuation is to shift the spinors of 
the negative helicity gluons 1 and 2 such that 

|T) = |l)+z|2), |2] = |2]-z|l]. (4.4) 

The corresponding momenta are also shifted, 

=P? + !<2|Y*|1], p£ -J#(*) =P£ " ^<2|7 M |l]. (4.5) 

We now consider the integral, 

1 f dz a , , 1 f dz , , , „ . s , , 
— d> —R n {z) = —& — (R n ( z ) - CR n {z)) . 4.6 
2m J c z 2m J c z 

Assuming that there is no surface term at infinity, the integral vanishes. The remaining 



residues are fixed by the multiparticle factorisation (4.1) so that the rational contribution 
is given by 

phys. poles z; 



^ A^\z)R R (z)+R L {z)Af{z) + ^ Reg ^ CR n (z) (4J) 



The last term is called the overlap term. It is quite simple to calculate if the poles at 
physical Zi are all first order, as this makes them similar to recursive terms. 

To sum up, the rational terms consist of the recursive terms which are similar in 
calculation to tree-level, completion terms which can be computed simply from the cut- 
containing functions, and the overlap terms which can be computed by considering the 
completion terms as certain internal momenta go on-shell. 
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(i + 1)- 




;+ 



i + (i + 1) 
2 - 1~ 





(i - 1)+ i+ 
1- T 

+ 





n 1 



(n-iy 




(i - ir 



Figure 6: The direct recursive diagrams contributing to R n (4>, 1 ,2 , . . . , n + ) with a |1|(2| shift. 



4.1 Recursive terms 



The recursive part of the rational contribution is denned by, 



D ^ Af{z)R R (z)+R L {z)Af(z) 

R n ~ 2^ ^2 



(4.8) 



For the choice of shift given in eq. (|4.4j ) the allowed types of contributing diagrams contri- 
butions are shown in fig. |6|. 



Combining the various diagrams, we find that recursive terms obey the following rela- 
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tion, 

R°(frl~,2-,3 + ,...,n + ) = 

71—1 1 

£ R&l-,^, (i + 1) + , • • • , n+)— AW(-P- p 2-, . . . , i+) 



i=4 



+ ^^°)(<A;T^P 2 7 J ,( i + l) + ,...,n + )—i?(-P+,2-, ...,.+) 

it — 1 _. 
+ ^i2(T-,-^+ ) i+,...,„+)_ A(°)(0;P < 7 1> 2-,...,(i-l)+) 

«=3 ' 
n— 1 1 

+ ^A(°)(1- -P f 7 lJ i+,...,n + )— i?(0;^+,2-,...,(z-l)+) 

i=4 Si ' 1 

+P(0; I" , P+ , 4+ , . . . , n+) — A<°> (-P^ , 2~ , 3+) 

+A(°)(T-, -P+,n + )-^P(<A; P n l,2- . . . , (n - 1) + ). (4.9) 

Here P represents the full rational part of the one- loop amplitude with fewer external legs. 
The contributions with i = n in the first term and i = 3 in the fourth term are absent 
because A2(</>, — , +) (and hence Pi2(</>> — , +)) is zero by conservation of angular momentum. 

Note that one could have written down other recursive terms that contain three-point 
pure gauge amplitudes. The most problematic of these are those where the "one- loop" - 
ness is in the gauge three-point, since the factorisation properties for complex momenta are 
not fully understood. However, when the two external gluons have opposite helicity it is 
guaranteed to vanish, because the corresponding splitting function does not have rational 
parts. For our particular choice of shift, the only possible pairings of the external gluons 
are (2~,3 + ) or (n + , 1~), and we find this significant simplification in all possible cases. 

This leaves the cases where pure gauge three-point is at tree level. Because of the way 
we have chosen our shift, the amplitudes 

A°(n+ I" -P- x ), A°(2- 3+ -P+ ) (4.10) 

both vanish, so that the only two contributions involving a three gluon vertex are when 
i = 3 in the first term and the last term of eq. (|4.9|). This latter contribution is the 
"homogenous" term in the recursion; it depends on the </>-MHV amplitude with one gluon 
fewer. The first two 0-MHV amplitudes are known, 

P 2 (0;l-,2-) = _La(°)(<M-,2-), (4.11) 

P 3 (0;l-,2-,3 + ) = -La( )(^,1-,2-,3+). (4-12) 

Because the tree amplitudes with fewer than two negative helicities vanish, the re- 
maining one-loop contributions needed are those with one negative helicity. These are 
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finite one-loop amplitudes and are entirely rational. The finite <j> — h • • • + amplitudes 
were computed for arbitrary numbers of positive helicity gluons in ref. [62]. As a concrete 
example, the three-gluon amplitude is given by, 

*<* = - ^4V;l-.2+,3+). (4.13) 

Similarly, the pure QCD — h • • • + amplitudes are given to all orders in ref. [49, 79]. In the 
four gluon case, the result is, 

ii 4 (l-,2+.3+,4+)^, , (4.14) 

4V ' ' ' ; 96vr 2 [12] (2 3) (34) [41] v ' 

The value that z takes is obtained by requiring that the shifted momenta 

^. = ^.±|(2| 7 ^|1], (4.15) 

is on-shell. In this equation, the sign is positive when the momentum set {pi,Pj} includes 
pi and is negative when it includes p 2 . There are four distinct channels, specified by the 
multiple invariants sn, s 2 , and the double invariants s n \ and s 2 3. In each channel, we find 
that the value of z and the hatted variables are given by, 



Si i channels: Z\ \ 



(2|Pi, 2 |l] 



|1> = -TW |2] = l^lp Pi - 1 = {2P,i\ ' (416) 



S2 i channels: z 2 j 



S2,i 



s n i channel: z n \ 



(2|P 3ji |l]' 

\iP3, 

(In) 



m _ 1-^1,^3,^2) |^i _ \P3,iPl,d] S _ |-Pl,il](2f3,i| (A-t 7 \ 

l ) ~ [lP^) ' |2J " (2P 2)i l] ' (2P 2>i l] ' ( j 



(2n)' 



|1} - |n) M' |2] "1^T' (2n) ' (418) 



S23 channel: z 23 



[32] 
[31] 



,T\_l P b23] |?1 _ hl [12] s |P 2 , 3 1][3| 

|1>_ W |2] - |3] pJ' 2 ' 3 "n^T" ( } 

4.2 Overlap terms 

The overlap terms are defined by [58] 

O n = J2^s z=Zi ^^-. (4.20) 

i 

They can be obtained by evaluating the residue of the cut completion term CR n given 
in eq. ( [3,28 ) in each of the physical channels. In order to make the residue calculation 



- 18 - 



straightforward, it is convenient to use identities such as — Sjj-i = {j\Pij-i\j] to 
rewrite eq. ( |3.28| ) in a way that exposes each of the physical poles, 

CRJ<f),l-,2-,3 + ,...,n + ) = 

V ' ' ' ' 12(nl)n™^(aa+l) 

(l|43|2>(l|A 3 4|2)f(l|A 3 4|2) - (14) [43] (3 2) 



S2 3 (4|P 2 , 3 |4] 2 




(l|nP 2 , n |2)(l|P 2 , n n|2)((l|P 2 , n n|2) 
+ v 


- (l|nP 2 , n |2)) 


S2,n(n\P2,n\n] 2 




(1|P 3 ,i3|2)(1|3P 3 ,i|2)((1|3P 3 ,i|2) - 

+ rr^ 


- <l]^3,i3|2>) 



(l|n(n - l)|2)(l|(n - l)P„,i|2) (l|(n - l)P n ,i|2) - {l\n(n - 1)|2) 



Sn i((n-l)|P n ,i|(n-l)] 2 



+ E 



i ! /(l|P i ,it|2)(l|iP i) i|2)((l|tP i)1 |2)-(l|Pi,it|2) 



(l|Pi,i(i - l)|2)(l|(i - l)P f ,i|2) ((l|(i - l)P,,i|2) - (l|P i(1 (i - 1)|2) 
+ - l)\P iA \(i ~ I)] 2 

! / r (l|tP2,«|2)(l|P2,ii|2)((l|P 2>i f|2)-(l|iP2,i|2) 

i=4 



S2,j \ S 2 ,i(i\P2,i\i 



2 



+ l)P 2ii |2)(l|P 2i i(z + 1)|2) (l|P2,i(i + 1)|2) - (l|(t + 1)P 2 , 4 |2) 



s 2j i((i + l)|P 2 ,i|(i + l)] 2 

(4-21) 

We observe that the cut completion term contains only simple residues so for the Pij 
pole, the overlap term is given by, 

OU = CR n ( Zij )^i (4.22) 

S i,j 

where Zij is the value of z that puts Pjj on-shell. The multiplicative factor removes the 
slj pole in CR n and replaces it with the correct propagator Sij. Note that the only terms 
that are affected by the momentum shifts are |1), |2] and any invariant including either pi 
or p2- The overall factor (nl) must be treated carefully, but not (23). 



Let us first consider the s 23 pole. The overlap term is given by, 
Ol 



-> 2 3 _ CR(Z 2 3)S23 



S23 

CrNp (14) [43] (T|P^4|2)((T|P^4|2) - (14) [43] (3 2) 
~ 2 [23] (nT)(4|P^ 3 |4]2n^- 2 (aa + l) 



(4.23) 
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Employing the definitions of the shifted variables in the s 23 channel given in eq. ( 4,19 ), we 
find that 



O 



2?> 



N P S i 23 [34](42)(4|P 12 |3]( Sl23 (42) + (4|P 12 |3]<32)) 



192^ S23 (n|P 12 |3](4|P 23 |l] 2 n™-i(aa + l) 
Similarly, we find that 

l2n N P s 1>n {2n)(2\P 3jn \n]{(2\P 3in P ltn \n) + s ljn (2n)) 



O 



O 



192^ S2i „(n|P ljn |l]2II^(aa + l) 

31 _ N P g3,n<32)(2|P 3 ,l|3](2|P3,li 3 3,n|3)((2|fh,ii^, n |3) + S 3 ,n<23)) 

~ 192vr 2 S3jl (n|P 3 , n P 3il |2)(3|P 3il |l] 2 n^- 2 (aa + l) 



(4.24) 

(4.25) 
(4.26) 



The overlap in the s n ± channel is complicated by the fact that there is an overall factor of 
l/(nl). We find, 



Ol 



n-l 

£ 

i=4 



(12)= 



192vr 2 (nl)Ulzl{aa + 1) 



[n(n-l)]((w-l)2) 
MP 



+ 



(n\P i>n (i - l)|2)(n|(i - l)P h2 \2)({n\(i - l)P t , 2 |2) - {n\P i>n (i - 1)|2)) 
(n|(i - l)P ij2 + Pi,„(i - 1)|2) 2 

1 1 



+ 



+ 



(n\P^ ljn P^ l>2 \2) (n\P itn Pi )2 \2) t 
(n|zP 3 ,,-i|2)(n|Pi,,_ii|2)((n|P M _ii|2) - (wji^j-ijg)) 

(n^i-l + P^-l^} 2 

1 1 



+ 



8,1 



(n|P l5i P 3ii |2) <n|Pi,i_ii%,i_i|2) 

Finally, when i < n — 1, the overlap terms are given by 
Np Si n 



(4.27) 



192tt 2 8 itl (n\P i>n P i>2 \2) KZl(aa + 1) 

(2|P i|2 (i - l)|2)(2|P i , 2 Pi,n|(i ~ l)){{2\Pi, 2 Pi,n\{i - 1)) + s hn (2(i - 1))) 

((i-l)\Pi,n\l? 
(2\P ia i\2)(2\P it2 P hn \i)((2\P h2 P hn \i) + s i|n <2i» 



+ 



O: 



2.1 



N P 



(4.28) 



1927T 2 S2i ,(n|P 1 ,,P 3 , i |2)n^- 2 (aa + l) 
(2t)(2|P3 i i|t](2|P3, i Pi, i |i)((2|P 3|i P 1 ,i|i) + ai,i(2i)) 



+ 



(^i,i|l] 2 

<2(t + l))(2|P 3ii |(i + l)]<2|P 3)i P 1 , i |(< + l))((2|P 3>i P lii |(i + 1)) + 8lji <2(< + 1))) 



{{i + i)\PxM 2 



(4.29) 
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5. The 4-point amplitude 



We will now compute the rational contribution explicitly for A^\(f> : 1~, 2~ , 3 + , 4 + ). The 
recursion relation (O) consists of four physical channels corresponding to poles in the 
invariants S23, S234, S41 and S341. 

Using the spinor shifts given in eq. (4.19) together with the known amplitudes given 
in eqs. fl4.11 )-( 4.14| ) it is straightforward to evaluate the direct recursive terms, 

i? 4 (^;l",2~,3 + ,4 + ) = P 234 + P 23 + Rf l + Rf . (5.1) 

For example, in the S234 channel, we find, 

Rf 4 = A (<M-,P 2 - 4 ) x — x i?(2-,3 + ,4 + ,-P+ 4 ) 

S234 



_ NpC T (1P 2 34) 2 (3P 23 4)[3P234] 3 
6 S234[23](34)(4P 234 )[P2342]' 

Employing the definitions of the shifts, ( [4.19| ), we easily find 



(5.2) 



D 234 _ 

N P [4 3]s 2 234 (2 4> 3 (3P 24 l] / - q\ 

4 967r 2 (34) 2 S2 34(2P34l](4P 2 3l] 2 ' 



Similarly, 



R 23 



N P s 123 [34][31](4|P 12 |3] 



4 ' 96vr 2 [23][12](4|P 23 |1] 2 8vr 2 



-^^°)(^,l-,2-,3 + ,4+), (5.4) 



i?! 1 =^ f A(°)(^l-,2-,3+,4+), (5.5) 
34! _ N P [31](2|P 13 |4] 3 

4 967r2 S341 (34)[41]2(2|P 3 4|l]- [ 1 

Notice that the expressions for P 23 and Pf 1 are not proportional to Np. It is clear that 
these two terms will cancel against similar terms appearing in the <^-MHV amplitude so 
that the Higgs amplitude does indeed have this property. 

The overlap terms are defined by, 

4 (</»;l-,2-,3 + ,4+) = of 4 + Of + Of 1 + Of . (5.7) 

The individual contributions can be obtained by setting n = 4 in eqs. ( [4.24 ), ( 4.26 ), (|4.2£| ) 
and ( [4.27 ) respectively. After some trivial simplifications, we find, 

23 N P Sl23 [34](42)( gl 23(42) + (4|P 12 |3](32)) 

4 ' 1927T 2 S2 3(34)(4|P 23 |1] 2 ' 1 j 

234 _ N P [34](42)m 2 J (m 2 / (42)-(2|P 34 Pi23|4)) 

4 1927T 2 S2 34(34)(4P 23 1] 2 ' 

n3 4i N P [34](2|P 13 |4]((2|P 13 |4] + (23)[34]) 

Ui ~ 192? s 34 i(34)[41] 2 ' (5 - iUj 

° A ~ 1927T 2 S4 i(34) • (5 - H) 
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With a little further algebra it is possible to write the sum of recursive and overlap 
terms in a form which is free of spurious singularities, 

£ 4 (0,l~,2-,3 + ,4+) = J R4(0;l",2-,3 + ,4+)+O 4 (</.;l-,2-,3+,4+) 



-L A W(A;l-,2-,3 + A + ) + 



iV P [43] 



8vr 2 v ' ' ' ' y 96tt 2 (34) 

(23)(1|P 24 |3] 2 | (41)(3|P 12 |3] (14)(2|P 13 |4] 2 | (32)(4|P 12 14] 



(34)[43][32]s 234 (34) [12] [32] (34)[43][41]s 3 4i (34) [12] [41] 
(12) 2 _ (12) _ (12)(2|P 13 |4] (12)(1|P 24 |3] (12) 2 (12) 2 1 
(34) [43] [12] 2[41]*34l 2[32]s 234 2s 23 2s 4 i J ' K ' 1 

where the pseudoscalar amplitude A^(A; 1~ ,2~ , 3 + , 4 + ) is given by the difference of 4> and 
cjr amplitudes. 

To form the full 1-loop <^>-MHV amplitude we need to add this to the cut-constructible 
piece given in eq. ( 3,24| ) and the completion term of eq. (4.21). 



4 1} ((M-,2-,3 + ,4 + ) = 
C 4 (0, 1", 2" 3 + , 4+)+CP 4 (0, 1", 2", 3+, 4+) + P 4 (0, 1", 2", 3 + , 4+). (5.13) 

As discussed earlier, the Higgs amplitude is constructed from the sum of the 4> and the 
$ amplitudes where the (jy contribution is obtained using parity symmetry, 

4 1) (F,l-,2-,3 + ,4+) = 4 1) (0,l-,2-,3 + ,4 + )+4 1) (0t,i~, 2 ~,3 + ,4 + ) 

= 4 1) (0,l^,2-,3 + ,4+) + (4 1) (0t,3-,4-,l+,2+)) 

(5.14) 

The one-loop amplitudes in this paper are computed in the four-dimensional helicity 
scheme and are not renormalised. To perform an MS renormalisation, one should subtract 
an MS counterterm from An , 

AM^AM-c r ^A£\ (5.15) 
The Wilson coefficient ( p.2j ) produces an additional finite contribution, 

4 1} -4 1} + y 4 0) - (5-16) 
6. Cross Checks and Limits 



6.1 Infra-red pole structure 

The infra-red poles are constrained to have a certain form proportional to the tree level 
amplitude [75, 76], 



n / 2 \ 6 
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Expanding the hyper geometric functions as a series in e quickly leads to a proof of this 
fact, and it can be seen that all logarithms vanish at order 1/e. 



6.2 Collinear Limits 

In general the collinear behaviour of one-loop amplitudes can be written[19, 80]: 

• • . » - lAl_1 ,P h ,i + 2 Al+2 , • • •) Splits (-P~ h ; i X \i + l Al+1 ) 

h=± 

• • ,i ~ 1 A<_1 ,P h ,i + 2 Al+2 , • • •) Split* 1 ) {-P~ h ; i x \i + l Al+1 ) (6.2) 



where the collinear limit is defined through pi — > zP and p^+i — > (1 — z)P. The universal 

splitting functions for QCD have been calculated in reference [19, 20, 81]. The tree- level 
splitting functions are given by, 

z 2 

Split (0) (-P + ,1~,2 + ) = -= (6.3) 

v^r^)<i2) 

Split* ) (-P+ 1+2") = ^£=£ (6.4) 





-z)(12) 


(1 




y/z(l 


-z){\2) 




1 



SpUt(°)(-P-,l+,2+) = (6.5) 

vn 1 " 2 )( 12 > 

Split*°)(-P-,r-,2-) = 0. (6.6) 

It is convenient to divide the one-loop splitting functions into cut-constructible and rational 
components, 

Split* 1 ) (-P- h , l Al , 2 Aa ) = Split* 1 )' C (-P- ft , l Al , 2 A2 ) + Split^'^-P^, l Al , 2 Aa ) (6.7) 
where 

Spht( 1 )' c (-P ± ,l-,2+) = Split(°)(-P ± ,l-,2+)Jx 

" ' f l- 2F 1 (l,-e;l-e;^-\- 2 F 1 (l,-e;l-e;—)), (6.8) 



S12/V V V z 



Split* 1 )' C (-P + ,l-,2-) = Split*°)(-P + ,l-,2-)^x 

,2 \ e 



/r S 1 - 2P1 (1, -e; 1 - e; ) - 2 P X ( 1, -e; 1 - e; — ) ) , (6.9) 



-512/ V V 

Split (1) ' c (-P~, I", 2") = 0, (6.10) 
Split* 1 )' R (-P ± , I", 2+) = 0, (6.11) 



Spht(^(-P + ,l-,2- ) = 7 ^^ 



96tt 2 [12] 



bplit t ^ ,1 ,2 j- g6vr2 [12]2 



(6.12) 
(6.13) 
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U - 1) 



(i + l) 



(i + 2) 
(j - 1) 




(i + 2) 

u - 1) 



Figure 7: Factorisation of combinations of two-mass "easy" box functions in limit where + 1. 



6.3 Collinear factorisation of the cut-constructible contributions 

Considering only the cut constructible contributions, we expect the following limiting be- 
haviour 

^(...Ai + i^,..)^ 1 

cfi (. . . , % - l^- 1 , P h , i + 2 A *+ 2 , . . .) Spilth (-P~ h ; i A * , i + l Al+1 ) 



h=± 



+A i ° ) _ 1 (...,i-l A *-\P h ,i + 2< 



i + 2 



-Split (1) ' c (-P- /l ;i A %i + l Al+1 ) (6.14) 



There are three collinear limits to consider for the 0-MHV amplitude; two collinear negative 
helicity gluons, two collinear positive helicity gluons and the mixed case with one of each 
helicity. 



We first note that the box and triangle scalar integrals always appear in the combina- 



tion 



i=l 

j n n+i— 1 ^ ,v 

- 2 X/ -^4 m6 ( s iJ) s i+lJ-i; s i,i+l) s i+l,i) ~ o X^ ^4 m ( g »,»+2; > s i+l,i 

i=l j=i+3 

and multiply the tree amplitude. 



+2j 



i=l 



(6.15) 



Many of the collinear factorisation properties of the individual box and triangle func- 
tions have been studied in Ref. [19] for example, 

F^ me (sij, Sj+ij_i; Sij-l, Sj+lj) + F^ me (si+ij, Si+2,j—l> s i+lJ-lj s i+2,j) 



i\\i+l „2me/ ^ 
— > 1*4 (Spj, Si+2,j-l'i s i+2,j, s P,j-l) 



(6.16) 



where spj = (P + Pi+2 + • • • + Pj) 2 an d which is illustrated in Fig. |7[ A similar relation 
applies when j = i + 3 such that the second and third terms are one-mass boxes. 

In much the same way, the particular combination of box and one-mass triangle func- 
tions shown in Fig. ||] factorise onto triangle functions as follows, 

F 2me (si, n , Si + 2,j-i; Si + 2,i) + nF 3 m (si jn ) — F 3 m (si4-i ) j_l) — F 3 lm (sj + 2,i) 



V (n-l)F^(s hn )-F^(s l+ 2,-i). 



(6.17) 
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(i + 1) 




Figure 8: Factorisation of a combination of two-mass "easy" box and one-mass triangle functions 
in limit where i\\i + 1. Note that the scale relevant for the triangle functions is indicated on the 
right of the triangle. The massless legs on the left of the triangles do not correspond to physical 
momenta. 

Taken together, it is straightforward to find the collinear limit of U n 
U n (l,...,i,i + l,...,n) -> U n -i(l,...,P,...,n) 

(6.18) 

which is independent of which (adjacent) pair are taken collinear. 

The remaining limits depend on the helicity of the collinear gluons. To explore the 
limits in detail, we find it convenient to utilise the expression for C n given in eq. ( |3.22| ). 

6.3.1 The mixed helicity collinear limit 

Let us now consider the case where we take a negative and a positive helicity collinear, 
e.g. the 2 1 1 3 limit shown in fig. |9|. We immediately see that the third diagram involves a 
factorisation onto a finite <f> amplitude which has no cut-constructible part and therefore 
gives no contribution. 

The finite logarithms also factorise in the expected way since, 

tr_(LP^ w (i - 1)2) 2 jj3 tr-(lP i|W (i - 1)P) 
su sip 
tr_(2P 3|i -i7l) 2j^3 tr_(PP 4 ,i-i»l) 
S12 sip 



i > 4, (6.19) 
i > 4, (6.20) 
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Figure 9: Collinear factorisation of A^\(f)\ 1,2, 3 + , . . . , n + ) taking p2 and p$ parallel, 
while 

tr_(lP 4 , n 32)) 2 ||3 



Sl2 

tr_(2341)) 2||3 

S12 



(6.22) 



ensures that terms with divergent logarithms, i.e. Lfc(s234, S23), will always be proportional 
to a trace which vanishes in the limit and hence do not appear in the one-loop splitting 
function. Together with eq. (6.18), we find that, as expected, 



C„(^;l-,2-,3 + ,...,n+) 2 ii 3 



A^ifr I", P", 4+, . . . , n+) Split«> C (-P+ 2~ 3+) 
+C n _!(</>; I" P- 4+, . . . , n+) Split(°)(-P+, 2" 3+). (6.23) 



6.3.2 Two collinear negative helicity gluons 

The 1||2 limit is rather trivial since there is only a single, finite term in this case as shown 
in fig. 10. The tree amplitude in this case vanishes, 

4 0) (l-,2-,3+,...,n+)^ 2 0, (6.24) 

therefore it is obvious that all the box and triangle terms of eq. fl3.24| ) will vanish. The 
remaining finite logs appear to have a singularity in S12, the worst coming from the traces 
raised to the 3rd power, 

tr_(LXY2) 3 _ (1|XY|2) 3 tr_(2XYl) 3 _ (2|XY|1) 3 

W2 " <12> 3 ' W2 " PP - ' (6 " 25) 

However the tree amplitude is proportional to (12) 3 so C n ((p, 1~, 2 _ , 3 + , . . . , n + ) vanishes 
in the limit as expected. 
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Figure 10: Collinear factorisation of (<p\ 1 ,2 ,3 + ,...,n + ) taking pi and P2 parallel. The 
only contribution is finite so the cut-constructible parts vanish in this limit. 

6.3.3 Two collinear positive helicity gluons 

The final limit occurs when we take any adjacent pair of positive helicities collinear. The 



relevant diagrams are shown in fig. 11. We can drop the third contribution since it involves 



the purely rational splitting function. 

We must also be able to show that there are no divergent terms coming from the finite 
logs and that these terms correctly factorise onto the lower point amplitude. This turns out 
to be slightly more involved than in the 2||3 case. First let us choose to take two adjacent 
particles a and b collinear where p a lies to the left of pt, in the clockwise ordering. Using, 

Sb,i a ~* (1 - z)sp,i + zs b+lji , (6.26) 

8 it a ZS i:P + (1 - Z)s it a-1, (6.27) 

it is then possible to show: 

tr_(2P 3 ,a-ial) fc r , > , tr_(2P 3 ,qfrl) fc T , , 

1 ^>k{ s 2,a, 82,a-l ) H J, L h\ s 2,b, S2,a) 

S 12 S 12 

a\\b tr_(2P 3a _iPl) fc . . 

% L k {S2,P,S2,a-l), (6.28) 

S 12 

and, 

tr_(lP 6+1 ,,62) fc tr_(lP M a2)\ 

J, L*k\ s b,ii s b+l,i) ' ^k{ s a,ij S 0i jJ 

s 12 s 12 

a\\b tr_(LP b+M P2) fc r 

-> 1 L k {sp ti ,s b+ i t i). (6.29) 

S 12 

Note that S2,p = (P2+- • -+Pa-i+P) 2 and spj = {P+p b+ i + . . .+pi) 2 . Using these identities 
and recognising that tr_(la62) and together with eq. ( |6.18| ), it is straightforward to 



show that eq. ( |3,24 ) has the correct factorisation properties, 

c n (</> ; i~,2-,3 + ,...,a+,& + ,...,n+)l fe 

C n _i(0; r, 2", 3+ . . . , a - 1+, P+, b + 1+, . . . , n+) Split(°)(-P~ a+, b+) 
+A i ° ) _ 1 ((f>; 1", 2~, 3+, . . . , a - 1+ P+, b + 1+, . . . , n+) Split«' C (-P-, a+, b + ). (6.30) 




(n-1)" 




Figure 11: Collinear factorisation of A^((f>; 1,2, 3 + , . . . , n + ) taking any two positive helicities 
parallel. 



6.3.4 Collinear factorisation of the rational contributions 

In this section, we focus on the collinear limit of the rational part of the 4 gluon amplitude 
with two negative and two positive helicities. There are three independent collinear limits, 
1112, 2113 and 31 14. 



Let us first consider the case where 1 and 2 are parallel. From fig. 10 we see 



that there is only a single contributing term and it is fairly straightforward to show that 
eq. ( fTp satisfies, 



1112 1 



i? 4 (^l-,2-,3+,4 + ) ± ^- I i?3(0 t ; J P-,3 + ,4+)Spht( o )(-P + ,l-,2-). 



3.31) 



The 2~||3 + limit, shown in fig. |9], has only two contributing diagrams for the raitonal 
piece since the first contribution vanishes since the one-loop mixed helicity splitting function 
has no rational part. The rational part of the four gluon amplitude can easily be shown to 
satisfy, 



i; I - , T , 3+ , 4+ ) + CR A (<P; 1~ , 2~ , 3+ , 4+) 
+i? 3 (0 t ;l", J P + ,4 + )Split(°)(-P-,2",3+) 
+R 3 (4>; r, P~, 4+) Split (°>(-P+ 2-, 3+) 



2 3 



(6.32) 



The final collinear limit is given by taking 3 + and 4 + collinear. Here we again find the 
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expected behaviour, 



i? 4 (^i-,2-,3+,4 + ) 3 i 4 
R 3 (<t>; 1~,2-,P + ) Split^ (-P- , 3 + , 4+) 
+^ (0) (0; 1" , 2" , P+) Split^'^-P" , 3+ , 4+) 

+^ (0) (</»; 1", 2~, P~) Split«' ii (-P + , 3+, 4+). (6.33) 



6.4 Soft Higgs Limit 



For the case of a massless Higgs boson, we can consider the kinematic limit pjj — > 0. 
Because of the form of the HG fJiU G^ u interaction, the Higgs field behaves like a constant 
in this limit, so the Higgs-plus-n-gluon amplitudes should be related to pure gauge theory 
amplitudes. Low energy theorems relate the amplitudes with zero Higgs momentum to 
pure gauge theory amplitudes [82] : 

A® (!!,{&,)*}) PH ^° Cg§-AV\{ 9i , AJ), (6.34) 

where C is the effective coupling of Higgs field to the gluon fields. The n-gluon tree 
amplitude is proportional to g n ~ 2 (see equation (|2.9|) ) therefore, 

4°) (H, { 9i , A*}) PH -Z° (const.) x (n-2)A^({ 9i ,\ i }). (6.35) 
In terms of the <fi and eft components [32], 

A^(^{ 9l ,X t }) (const.) x (n_ - 1)A^ ({ 9t , X t }), (6.36) 
4 0) (^,{9i, M) P C° (const.) x ( n+ - l)A^({ 9i , A*}) (6.37) 

where n+ and are the number of positive and negative helicity particles respectively. 

The one-loop amplitudes are proportional to g n hence similarly one can deduce the 
following behaviour in the soft Higgs limit, 

A^{H,{ 9l ,\}) PH ^° (const.) x nA${{9i,\})- (6-38) 

It has been conjectured in [62] that the <f> and (jr components generalise from the tree level 
relations to give, 

4 1 )(0, { 9i , A,}) P C° (const.) x n-AV)({gi, AJ), (6.39) 
A^(^,{ 9u X t }) P *£° (const.) x n + A^ ({ 9l , X t }). (6.40) 
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The 4-gluon MHV amplitude at one-loop in QCD has been derived in ref. [83] and is 
given, unrenormalised in the four-dimensional helicity scheme by 



C 4 (l-, 2", 3+ 4+) = - ^#(r, 2" 3+ 4+) 



2 \ e 



-Si4 



+ 



/3b 



iVe(l - 2e) V-s 



2 \ e 



'14 



R41-, 2-, 3+ 4+) =^A(°) (I", 2T, 3+, * 



^Y + log 2 ^)^ 2 

-512/ V s 14, 



(6.41) 
(6.42) 



6.4.1 Soft limit of A 



(!)/ 



i,l-,2-,3+ 4 H 



Firstly let us consider the soft limit of the cut constructible components, eq. ( |3.24 ). The 
1-mass and 2-mass easy box functions and triangle functions have smooth soft limits where 
we take: 







fx 



2 \ 6 



0. 



(6.43) 
(6.44) 



We must apply the same relations to the finite logs coming from the tensor triangle integrals, 
for instance we find: 



£fc(s34l> s 4l) 



Bub(s 34 i) - Bub(s 4 i) P4 ,^o (-1) 



(•S341 — S41) 



sf A e(l - 2e) \-s u 



2 \ e 



(6.45) 



Applying these relations together with momentum conservation to the cut-constructible 
part of the 4 gluon amplitude we find that the boxes and triangle functions collapse onto the 
correct 1/e 2 poles while the coefficients of the bubble contributions simplify considerably, 
e.g., 



S12 J 



14' 



so that all of the logarithms proportional to Np cancel among themselves leaving, 

.2 \ e 



2/?o 



s M y iVe(l-2e) 
Combining this with the 1/e 2 poles we find, 

C4(0,l-,2-,3 + ,4 + ) P -^ O 2C4(l-,2-,3 + ,4 + ) 
The soft limit of the rational part, eq. G5.12|) and eq. ( |3.28| ), gives 

i?4(^;l-,2-,3+,4+) P n -^A(°)(l-,2-,3 + ,4+), 



(6.46) 



(6.47) 



(6.48) 



(6.49) 
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with all the poles in the triple invariants vanishing. This implies that the rational terms 
vanish in the limit for the full Higgs field, 

i? 4 (^;l _ ,2-,3 + ,4 + ) P ^°0, (6.50) 

since the $ contribution appears with the opposite sign. This confirms that cut-constructible 
terms of the amplitude do follow the naive factorisation proposed in [62] but the rational 
parts do not. This seems to be consistent with the results for the HqqQQ amplitude 
presented in [64]. 

7. Conclusions 

Recent developments based on unitarity and on-shell recursion relations have lead to impor- 
tant breakthroughs in computing non-super symmetric one-loop gauge theory amplitudes. 
In particular, new compact analytic results have been obtained for gluonic amplitudes 
involving six or more gluons. 

In this paper we have focused on amplitudes involving gluons and a colourless scalar 
- the Higgs boson. The model which we use to calculate these amplitudes is the tree-level 
pure gauge theory plus an effective interaction HG fJ , u G fll/ produced by considering the heavy 
top quark limit of the one-loop coupling of Higgs and gluons in the non-supersymmetric 
standard model. Following Ref. [32], we split the interaction into selfdual and anti-selfdual 
pieces. The self-dual (anti-self-dual) gauge fields interact with the 4> (0') scalars respec- 
tively and because of this selfduality, the amplitudes for (f> plus n gluons, and those for $ 
plus n gluons, each have a simpler structure than the gluonic amplitudes for either H or 
A. 

Previous studies using the unitarity-factorisation booststrap program have focused on 
the (f) amplitudes that are finite at one-loop, i.e. the all positive or nearly all positive 
helicity [62] or the divergent amplitudes when the gluons all have negative helicity [63]. In 
this paper we have employed four-dimensional unitarity and recursion relations to compute 
the one-loop corrections to amplitudes involving a colourless scalar <f>, two colour adjacent 
negative helicity gluons and an arbitrary number of positive helicity gluons - the so-called 
0-MHV amplitudes. 

The gluonic production of Higgs bosons via a heavy quark loop is expected to be the 
largest source of Higgs bosons at the LHC. Because of the size of the strong coupling, it will 
be important to understand Higgs plus jet events in some detail. The two jet channel has 
contributions of both QCD [72, 73, 84-89] and weak boson fusion origin [90-92]. Separating 
these two "signals" is of crucial importance for measuring the coupling of the Higgs bosons 
to standard model particles from LHC data [93]. The amplitudes presented here may be 
useful in computing the gluon fusion contamination of the weak boson fusion signal. 

The one-loop amplitudes naturally divide into cut-containing C n and rational parts 
R n . We used the double cut unitarity approach of ref. [37] to derive all the multiplicity 
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results for C n given in eq. ( |3.24 ). The rational terms have several sources - first the cut- 



completion term CR n which eliminates the unphysical poles present in C n , second the 
direct on-shell recursion contribution and third the overlap term O n . Explicit formulae 
for these contributions are given in eqs. ( |3.28| ), ( f4.9| ) and (4.24)— ( 4.29| ). An explicit solution 



for the four gluon case is given in section [|. We have checked our results in the limit where 
two of the gluons are collinear, in the limit where the scalar becomes soft and against 
previously known results for up to four gluons. 
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A. Spinor conventions 

In the spinor helicity formalism [94-99] an on-shell momentum of a massless particle, 
k^k 11 = 0, is represented as 

where X a and A^ are two commuting spinors of positive and negative chirality. Spinor 
inner products are defined by 

(A, A') = e a p\ a X'^ , [A, A'] = -e^X'* , (A.2) 

and a scalar product of two null vectors, k aa = X a Xa and p a a = A' Q A^, becomes 

k^= -i(A,A')[A,A']. (A.3) 

We use the shorthand (ij) and [ij] for the inner products of the spinors corresponding to 
momenta pi and pj , 

(ij) = (X l ,X j ), [ij] = [X i ,X j ]. (A A) 
For gluon polarization vectors we use 

ej(k,Q = ± { ^ k l ) , (A.5) 



where k is the gluon momentum and £ is the reference momentum, an arbitrary null vector 
which can be represented as the product of two reference spinors, ^ a a = CoCa- 



- 32 - 



B. Scalar integrals 



The one-loop functions that appear in the all-orders cut-constructible contribution C n given 
in eq. ( [3.24 ) are defined by, 



F 2me (P 2 ,Q 2 ;s,t) = ^ 



+ 



Bub(s) 



2 \ e 



2 F X fl,- e ;l-e;-| 



+ ( r_ ] 2 F X ( l,- e; l- e; -- 



4) ^fl,-,;!-,;-! 



2 \ e 



2 \ e 



F X l,-e;l-e; 



_p2 

^ ! \--<:\-r. 



2 \ e 
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~st 

US 



P 2 Q 2 - st 
ut 



2 \ e 



p2 

2 \ e 



Fi l,-e;l-e; 



uP 2 



Q 2 

1 



2 Fi l,-e;l-e; 



F 2 Q 2 - st 
F 2 Q 2 - 



e(l - 2e) V-8 



2 \ e 



(B.l) 



(B.2) 



(B.3) 
(B.4) 
(B.5) 



References 

[1] M. L. Mangano, M. Moretti and R. Pittau, Multijet matrix elements and shower 
evolution in hadronic collisions: W b anti-b + (n)jets as a case study, Nucl. Phys. 
B632 (2002) 343-362 [|hep-ph/0108069|| . 



[2] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, Alpgen, a 
generator for hard multiparton processes in hadronic collisions, JEEP 07 (2003) 001 
flhep-ph/0206293[| . 



[3] A. Pukhov et. al, Comphep: A package for evaluation of feynman diagrams and 
integration over multi-particle phase space, user's manual for version 33, 
tiep-ph/9908288; . 



- 33 - 



CompHEP Collaboration, E. Boos et. al., Comphep 4-4 : Automatic computations 
from lagrangians to events, Nucl. Instrum. Meth. A534 (2004) 250-259 
[frep-ph/0403113fl . 

A. Kanaki and C. G. Papadopoulos, Helac: A package to compute electroweak 



helicity amplitudes, Comput. Phys. Commun. 132 (2000) 306-315 [frep-ph/0002082(1 . 



C. G. Papadopoulos, Phegas: A phase space generator for automatic cross-section 
computation, Comput. Phys. Commun. 137 (2001) 247-254 [|hep-ph/0007335|1 . 



T. Stelzer and W. F. Long, Automatic generation of tree level helicity amplitudes, 



Comput. Phys. Commun. 81 (1994) 357-371 lhep-ph/9401258 ]. 



F. Maltoni and T. Stelzer, Madevent: Automatic event generation with madgraph, 



JEEP 02 (2003) 027 jhep-ph/0208156 |. 



T. Gleisberg et. al., Sherpa 1. alpha, a proof- of- concept version, JEEP 02 (2004) 056 
[frep-ph/031 1263fl . 

F. Krauss, R. Kuhn and G. Soft", Amegic++ 1.0: A matrix element generator in 
C++, JEEP 02 (2002) 044 fcep-ph/0 109036(1 . 



F. Caravaglios, M. L. Mangano, M. Moretti and R. Pittau, A new approach to 
multi-jet calculations in hadron collisions, Nucl. Phys. B539 (1999) 215-232 
jhep-ph/9807570| |. 

F. A. Berends and W. T. Giele, Recursive calculations for processes with n gluons, 
Nucl. Phys. B306 (1988) 759. 

D. A. Kosower, Light cone recurrence relations for QCD amplitudes, Nucl. Phys. 
B335 (1990) 23. 

P. Draggiotis, R. H. P. Kleiss and C. G. Papadopoulos, On the computation of 



multigluon amplitudes, Phys. Lett. B439 (1998) 157-164 ltiep-ph/9807207 



S. Catani and M. H. Seymour, A general algorithm for calculating jet cross sections 
in NLO QCD, Nucl. Phys. B485 (1997) 291-419 jhep-ph/9605323 l. 



R. K. Ellis, W. T. Giele and G. Zanderighi, The one-loop amplitude for six-gluon 
scattering, JEEP 05 (2006) 027 1 |hep-ph/0602185| ] . 



A. Denner, S. Dittmaier, M. Roth and L. H. Wieders, Electroweak corrections to 
charged- current e + e~ — ► 4 fermion processes: Technical details and further results, 
Nucl. Phys. B724 (2005) 247-294 |frep-ph/0505042|l . 



A. Denner, S. Dittmaier, M. Roth and L. H. Wieders, Complete electroweak \{a) 
corrections to charged- current e + e~ — > 4 fermion processes, Phys. Lett. B612 (2005) 
223-232 1frep-ph/0502063|l . 



- 34 - 



Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, One loop n-point gauge 
theory amplitudes, unitarity and collinear limits, Nucl. Phys. B425 (1994) 217-260 
[|hep-ph/9403226(| . 



Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Fusing gauge theory tree 



amplitudes into loop amplitudes, Nucl. Phys. B435 (1995) 59-101 [ tiep-ph/9409265 1 



E. Witten, Perturbative gauge theory as a string theory in twistor space, Commun. 
Math. Phys. 252 (2004) 189-258 |hep-th/0312171||. 



F. Cachazo, P. Svrcek and E. Witten, MEV vertices and tree amplitudes in gauge 



theory, JEEP 09 (2004) 006 flhep-th/0403047 . 

R. Britto, F. Cachazo and B. Feng, New recursion relations for tree amplitudes of 



gluons, Nucl. Phys. B715 (2005) 499-522 ||hep-th/0412308 |. 

R. Britto, F. Cachazo, B. Feng and E. Witten, Direct proof of tree-level recursion 



relation in Yang-Mills theory, Phys. Rev. Lett. 94 (2005) 181602 |hep-th/0501052 |. 



K. Risager, A direct proof of the CSW rules, JEEP 12 (2005) 003 [|hep-th/0508206|1 . 



S. D. Badger, E. W. N. Glover, V. V. Khoze and P. Svrcek, Recursion relations for 
gauge theory amplitudes with massive particles, JEEP 07 (2005) 025 
[|hep-th/0504159|| . 



D. Forde and D. A. Kosower, All-multiplicity amplitudes with massive scalars, Phys. 
Rev. D73 (2006) 065007 [[hep-th/0507292l. 



S. D. Badger, E. W. N. Glover and V. V. Khoze, Recursion relations for gauge 
theory amplitudes with massive vector bosons and fermions, JEEP 01 (2006) 066 



[|hep-th/050716l l. 

C. Schwinn and S. Weinzierl, On-shell recursion relations for all born QCD 



amplitudes, [hep-ph/ 0703021 



C. Schwinn and S. Weinzierl, SUSY ward identities for multi-gluon helicity 



amplitudes with massive quarks, JEEP 03 (2006) 030 |hep-th/0602012 



P. Ferrario, G. Rodrigo and P. Talavera, Compact multigluonic scattering amplitudes 
with heavy scalars and fermions, Phys. Rev. Lett. 96 (2006) 182001 
[|hep-th/0602043(| . 



L. J. Dixon, E. W. N. Glover and V. V. Khoze, MEV rules for Eiggs plus 
multi-gluon amplitudes, JEEP 12 (2004) 015 [hep-th/0411092 ]. 



S. D. Badger, E. W. N. Glover and V. V. Khoze, MEV rules for Eiggs plus 



multi-parton amplitudes, JEEP 03 (2005) 023 jhep-th/0412275 l 



- 35 - 



Z. Bern, D. Forde, D. A. Kosower and P. Mastrolia, Twistor-inspired construction of 



electroweak vector boson currents, Phys. Rev. D72 (2005) 025006 [ tiep-ph/0412167 1. 



F. Cachazo and P. Svrcek, Lectures on twistor strings and perturbative Yang-Mills 



theory, PoS RTN2005 (2005) 004 flriep-th/0504194|1 . 



Z. Bern, L. J. Dixon and D. A. Kosower, On-shell methods in perturbative QCD, 
|arXiv: 0704 .2798 [hep-ph]| . 

A. Brandhuber, B. Spence and G. Travaglini, One-loop gauge theory amplitudes in 
M = 4 super Yang-Mills from MHV vertices, Nucl. Phys. B706 (2005) 150-180 



tiep-th/0407214 



R. Britto, F. Cachazo and B. Feng, Generalized unitarity and one-loop amplitudes in 
M = 4 super Yang-Mills, Nucl. Phys. B725 (2005) 275-305 [|hep-th/0412103|l . 



R. Britto, F. Cachazo and B. Feng, Computing one-loop amplitudes from the 
holomorphic anomaly of unitarity cuts, Phys. Rev. D71 (2005) 025012 
Hriep-tri/0410179|| . 



R. Britto, E. Buchbinder, F. Cachazo and B. Feng, One-loop amplitudes of gluons in 
SQCD, Phys. Rev. D72 (2005) 065012 ^ep-ph/050"313^| . 



R. Britto, B. Feng and P. Mastrolia, The cut-constructible part of QCD amplitudes, 



Phys. Rev. D73 (2006) 105004 [|hep-ph/0662178 l 



J. Bedford, A. Brandhuber, B. Spence and G. Travaglini, N on- super symmetric loop 
amplitudes and MHV vertices, Nucl. Phys. B712 (2005) 59-85 [|hep-th/0412108|1 . 



D. Forde, Direct extraction of one-loop integral coefficients, arXiv : 0704 . 1835 



[hep-ph] 



W. L. van Neerven, Dimensional regularization of mass and infrared singularities in 
two loop on-shell vertex functions, Nucl. Phys. B268 (1986) 453. 

C. Anastasiou, R. Britto, B. Feng, Z. Kunszt and P. Mastrolia, D-dimensional 



unitarity cut method, Phys. Lett. B645 (2007) 213-216 |hep-ph/060919l |. 



P. Mastrolia, On triple-cut of scattering amplitudes, Phys. Lett. B644 (2007) 272 
[arXiv:hep-th/0611091]. 

R. Britto and B. Feng, Unitarity cuts with massive propagators and algebraic 



expressions for coefficients, tiep-ph/0612089. 



C. Anastasiou, R. Britto, B. Feng, Z. Kunszt and P. Mastrolia, Unitarity cuts and 
reduction to master integrals in d dimensions for one-loop amplitudes, JHEP 03 
(2007) 111 [|hep-ph/0612277]| . 



- 36 - 



Z. Bern, L. J. Dixon and D. A. Kosower, On-shell recurrence relations for one-loop 
QCD amplitudes, Phys. Rev. D71 (2005) 105013 friep-th/0501240 |. 



Z. Bern, L. J. Dixon and D. A. Kosower, The last of the finite loop amplitudes in 



QCD, Phys. Rev. D72 (2005) 125003 lhep-ph/0505055 1 



Z. Bern, L. J. Dixon and D. A. Kosower, Bootstrapping multi-parton loop amplitudes 
QCD, Phys. Rev. D73 (2006) 065013 [ |hep-ph/0507005| ]. 



in 



Z. Xiao, G. Yang and C.-J. Zhu, The rational part of QCD amplitude. I: The general 
formalism, Nucl. Phys. B758 (2006) 1-34 |hep-ph/0607015|| . 



Z. Xiao, G. Yang and C.-J. Zhu, The rational part of QCD amplitude. Ill: The 
-gluon, Nucl. Phys. B758 (2006) 53-89 lhep-ph/0607017|| . 



six- 



G. Ossola, C. G. Papadopoulos and R. Pittau, Reducing full one-loop amplitudes to 
scalar integrals at the integrand level, Nucl. Phys. B763 (2007) 147 
[arXiv:hep-ph/0609007]. 

T. Binoth, J. P. Guillet and G. Heinrich, Algebraic evaluation of rational polynomials 
in one-loop amplitudes, JEEP 02 (2007) 013 ]hep-ph/0609054 l. 



X. Su, Z. Xiao, G. Yang and C.-J. Zhu, The rational part of QCD amplitude. E: The 
five-gluon, Nucl. Phys. B758 (2006) 35-52 |hep-ph/0607016|1 . 



D. Forde and D. A. Kosower, All-multiplicity one-loop corrections to MEV 



amplitudes in QCD, Phys. Rev. D73 (2006) 061701 [hep-ph/0509358]. 



C. F. Berger, Z. Bern, L. J. Dixon, D. Forde and D. A. Kosower, Bootstrapping 
one-loop QCD amplitudes with general helicities, Phys. Rev. D74 (2006) 036009 
[|hep-ph/0604195(| . 

C. F. Berger, Z. Bern, L. J. Dixon, D. Forde and D. A. Kosower, All one-loop 
maximally helicity violating gluonic amplitudes in QCD, Phys. Rev. D75 (2007) 



016006 flhep-ph/0607014 



T. Binoth, G. Heinrich, T. Gehrmann and P. Mastrolia, Six-Photon Amplitudes, 
arXiv:hep-ph/0703311. 

G. Ossola, C. G. Papadopoulos and R. Pittau, Numerical Evaluation of Six-Photon 
Amplitudes, arXiv:0704.1271 [hep-ph]. 

C. F. Berger, V. Del Duca and L. J. Dixon, Recursive construction of 
higgs+multiparton loop amplitudes: The last of the (fi-nite loop amplitudes, Phys. 



Rev. D74 (2006) 094021 [|hep-ph/0608180 



S. D. Badger and E. W. N. Glover, One-loop helicity amplitudes for E — > gluons: 
The all- minus configuration, Nucl. Phys. Proc. Suppl. 160 (2006) 71-75 
[|liep-ph/0607 139fl . 



-37- 



[64] R. K. Ellis, W. T. Giele and G. Zanderighi, Virtual QCD corrections to higgs boson 
plus four parton processes, Phys. Rev. D72 (2005) 054018 |frep-ph/0506196|l . 



[65] M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and V. I. Zakharov, Low-energy 
theorems for Higgs boson couplings to photons, Sou. J. Nucl. Phys. 30 (1979) 
711-716. 

[66] F. Wilczek, Decays of heavy vector mesons into Higgs particles, Phys. Rev. Lett. 39 
(1977) 1304. 

[67] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, Decoupling relations to C(af ) 
and their connection to low-energy theorems, Nucl. Phys. B510 (1998) 61-87 
[|hep-ph/9708255[] . 



[68] T. Inami, T. Kubota and Y. Okada, Effective gauge theory and the effect of heavy 
quarks in Higgs boson decays, Z. Phys. C18 (1983) 69. 

[69] U. Baur and E. W. N. Glover, Higgs boson production at large transverse momentum 
in hadronic collisions, Nucl. Phys. B339 (1990) 38-66. 

[70] R. K. Ellis, I. Hinchliffe, M. Soldate and J. J. van der Bij, Higgs decay to t + t~ : A 
possible signature of intermediate mass higgs bosons at the ssc, Nucl. Phys. B297 
(1988) 221. 

[71] M. Kramer, E. Laenen and M. Spira, Soft gluon radiation in higgs boson production 



at the LHC, Nucl. Phys. B511 (1998) 523-549 jhep-ph/96 11272 



[72] V. Del Duca, A. Frizzo and F. Maltoni, Higgs boson production in association with 



three jets, JHEP 05 (2004) 064 flhep-ph/0404013 



[73] S. Dawson and R. P. Kauffman, Higgs boson plus multi - jet rates at the SSC, Phys. 
Rev. Lett. 68 (1992) 2273-2276. 

[74] S. J. Parke and T. R. Taylor, An amplitude for n gluon scattering, Phys. Rev. Lett. 
56 (1986) 2459. 

[75] S. Catani, The singular behaviour of QCD amplitudes at two-loop order, Phys. Lett. 
B427 (1998) 161-171 [[hep-ph/9802439" |. 



[76] W. T. Giele and E. W. N. Glover, Higher order corrections to jet cross-sections in 
e + e~ annihilation, Phys. Rev. D46 (1992) 1980-2010. 

[77] C. R. Schmidt, h — > ggg(gqq) at two loops in the large-m(t) limit, Phys. Lett. B413 
(1997) 391-395 [ frep-ph/9707448]] . 

[78] Z. Bern and G. Chalmers, Factorization in one loop gauge theory, Nucl. Phys. B447 
(1995) 465-518 [frep-ph/950323611 . 



- 38 - 



G. Mahlon, Multi - gluon helicity amplitudes involving a quark loop, Phys. Rev. D49 



(1994) 4438-4453 jhep-ph/9312276 



D. A. Kosower, All-order collinear behavior in gauge theories, Nucl. Phys. B552 
(1999) 319-336 [tiep-ph/990120l[| . 



Z. Bern, V. Del Duca, W. B. Kilgore and C. R. Schmidt, The infrared behavior of 
one-loop QCD amplitudes at next-to-next-to-leading order, Phys. Rev. D60 (1999) 
116001 [|hep-ph/9903516 [. 



J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter's Guide. 
Addison-Wesley, 1990. 

Z. Bern and D. A. Kosower, The computation of loop amplitudes in gauge theories, 
Nucl. Phys. B379 (1992) 451-561. 

V. Del Duca, W. Kilgore, C. Oleari, C. R. Schmidt and D. Zeppenfeld, Kinematical 
limits on Higgs boson production via gluon fusion in association with jets, Phys. Rev. 
D67 (2003) 073003 [|hep-ph/030lol3]1 . 



R. P. Kauffman, S. V. Desai and D. Risal, Production of a Higgs boson plus two jets 
in hadronic collisions, Phys. Rev. D55 (1997) 4005-4015 |iep-ph/961054ll . 



V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and D. Zeppenfeld, H + 2 jets via 



gluon fusion, Phys. Rev. Lett. 87 (2001) 122001 [tiep-ph/0105129 



V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and D. Zeppenfeld, Gluon-fusion 
contributions to H + 2 jet production, Nucl. Phys. B616 (2001) 367-399 
|hep-pri/0108030(| . 



G. Klamke and D. Zeppenfeld, Higgs plus two jet production via gluon fusion as a 



signal at the CERN LHC, |hep-pri/0703202 . 



J. M. Campbell, R. Keith Ellis and G. Zanderighi, Next-to-leading order higgs + 2 
jet production via gluon fusion, JHEP 10 (2006) 028 |hep-ph/0608194]. 



E. L. Berger and J. Campbell, Higgs boson production in weak boson fusion at 



next-to- leading order, Phys. Rev. D70 (2004) 073011 jhep-ph/0403194 



T. Figy, C. Oleari and D. Zeppenfeld, Next-to-leading order jet distributions for higgs 
boson production via weak-boson fusion, Phys. Rev. D68 (2003) 073005 
|hep-pri/0306109|| . 

T. Figy and D. Zeppenfeld, QCD corrections to jet correlations in weak boson fusion, 
Phys. Lett. B591 (2004) 297-303 |hep-ph/0403297 l. 



M. Duhrssen et. at, Extracting higgs boson couplings from Ihc data, Phys. Rev. D70 



(2004) 113009 flhep-ph/0406323 1 



- 39 - 



[94] F. A. Berends, R. Kleiss, P. De Causmaecker, R. Gastmans and T. T. Wu, Single 
bremsstrahlung processes in gauge theories, Phys. Lett. B103 (1981) 124. 

[95] P. De Causmaecker, R. Gastmans, W. Troost and T. T. Wu, Multiple bremsstrahlung 
in gauge theories at high- energies. 1. general formalism for quantum electrodynamics, 
Nucl. Phys. B206 (1982) 53. 

[96] R. Kleiss and W. J. Stirling, Spinor techniques for calculating p anti-p — > w^/z + 
jets, Nucl. Phys. B262 (1985) 235-262. 

[97] J. F. Gunion and Z. Kunszt, Improved analytic techniques for tree graph calculations 
and the g g q anti-q lepton anti-lepton subprocess, Phys. Lett. B161 (1985) 333. 

[98] Z. Xu, D.-H. Zhang and L. Chang, Helicity amplitudes for multiple bremsstrahlung in 
massless nonabelian gauge theories, Nucl. Phys. B291 (1987) 392. 

[99] M. L. Mangano and S. J. Parke, Multiparton amplitudes in gauge theories, Phys. 
Rept. 200 (1991) 301-367. 



-40- 



